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Abstract 

We derive simple expressions that relate the noise and correlation properties of a general time-dependent quantum 
conductor to the wave functions of the system. The formalism provides a practical route for numerical calculations 
of quantum noise in an externally driven system. We illustrate the approach with numerical calculations of the noise 
properties associated to a voltage pulse applied on a one-dimensional conductor. The methodology is hower fully 
general and can be used for a large class of mesoscopic conductors. 


1. Introduction 

Among the many contributions of Markus Biittiker to the field of mesoscopic physics (now best known as nano¬ 
electronics), his pioneer work on time-dependent phenomena was particularly dear to him. He was insistent in pointing 
out the role of displacement currents (required to restore current conservation), the need for a theory that preserved 
“gauge invariance” and more generally built up the general framework and concepts to address this physics. While 
the key theoretical works were performed in the early 90s III] mail a, the corresponding experiments were diffi¬ 
cult (GHz physics at mK temperatures) so that more than 10 years elapsed before the hrst quantum RC circuit could 
actually be measured El. The held has considerably matured since then, with the latest generation of experiments 
performed directly in the time domain Q. 

Markus visited Saclay on March 2006 which is when the senior author of this paper was introduced to the subject. 
At that time, although the analytical theory was well developed, its computational counterpart was still in its infancy 
(see m for a short history). We now have very effective numerical tools ISl, with a computational effort linear both 
in time and system size, so that numerical complexity is no longer an issue 0 to simulate time-dependent systems. 
In this article, we would like to extend these tools to calculate another quantity, also very dear to Markus, namely the 
quantum fluctuations of observables. Quantum noise is a quantity which is not only sensitive to the wave aspect of 
quantum transport but also to particle statistics QO). We shall see that its numerical calculation requires some care in 
order to disentangle the relevant contributions from the large background already present at equilibrium. 

2. The Pauli principle in driven mesoscopic systems 

Let us consider a general mesoscopic system described by a time dependent (externally driven) quadratic Hamil¬ 
tonian, 

( 1 ) 

Uj 

where cj (cj) are the usual Fermionic creation (annihilation) operators of a one-particle state on site i. The site index 
i typically labels space as well as spin, orbital (s,p,d,f) and/or Nambu (for superconductors) degrees of freedom. The 
system is open, i.e. consists of a finite time-dependent central part connected to infinite (stationary) electrodes. We 
further suppose that the Hamiltonian is time independent for f < 0 and we switch on the time-dependent part (voltage 
pulses, light...) at f > 0. 
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For f < 0, the solution to this problem is well known: one introduces the stationary scattering states (labeled 
by their energy E and lead mode a) which diagonalize the one-body Hamiltonian, 

^ EXUO (2) 

J 

and build many-body Slater determinants from these one-body states. For quantities such as the current, this amounts 
to filling up the energy E with probability fa{E) (Fermi function of the lead to which channel a belongs) which leads 
to the celebrated Landauer formula for the conductance. Within our notations, the current between sites i and j reads, 

a 

This approach is by now so standard that can be obtained directly from open source softwares such as the Kwant 
package im. The generalization of the above approach to time-dependent phenomena is in fact straightforward: one 
simply follows the evolution of the scattering states upon switching on the time-dependent perturbation and solves, 

ih-'i'aEii, f) = X Ho<f)'F«£(7, t) (4) 

j 

with the initial condition 

'F„£(/,f = 0) = 'F^'£(0. (5) 

The time-dependent current is now given by, 

hjit) = -y Im ^ ^ dE /„(£)['F„£(/, f)]*H,;(f)'F,£(J, t). (6) 

Equation Q is transparent physically but not directly useful for numerical computations as the wave-function '^aEiU 0 
is a vector of infinite size. However, by studying the deviation between 'i'^EiU 0 and 'F"g(/)e"'^' one obtains a finite 
vector amenable to a numerical solution jSl. In order to prove that Equations Q, Q and (|^ lead to the correct 
generalization of the Landauer formula. Ref. IS) proved its mathematical equivalence with the Wingreen-Meir ifT^ 
approach based on the Keldysh formalism. One can also prove that the 'i'aEiU 0 are well defined scattering states, i.e. 
have the correct structure of superposition of incoming and outgoing states in the leads. The Pauli principle is fully 
enforced within this scheme, as the initial orthogonality relations of the states are preserved by the unitary evolution of 
Schrodinger equation. It is interesting to contrast the above set of equations with studies where a (usually Gaussian) 
wavepacket is propagated through the system. In fine, the dynamical equation solved is actually the same. There 
are however two crucial differences: the initial boundary condition (delocalized overt the infinite system in our case) 
and the final integration over energy which restores the fermionic statistics. The approach presented here is fully 
many-body and treats Pauli principle exactly. 

3. Quantum Noise 

We now discuss how the approach outlined above can be generalized to calculate the noise properties of a general 
mesoscopic system subject to time-dependent perturbation. We will focus in particular on the variance of the total 
number of particles sent through an electrode fi. This quantity is of special interest to us as it is relatively easy to 
measure experimentally (upon sending repeated pulses it is essentially a d.c. measurement as opposed to much more 
challenging high frequency measurements) and it is conserved and gauge invariant in Markus’ sense |[8l. 

3.1. General expressions 

The current operator is defined as a sum of the local currents flowing through a cross section corresponding to an 
electrode p (in practice p is the collection of hopping elements that connect the central system to the electrode) 

4(f) = ^ EIp{t)c]{t)cj{t) - nji{t)c]{t)ci(t), (7) 
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and we define the current-current correlation function as 


S,y(t, 0 = {lit) - (lit))] X {lit') - (lit'))) . (8) 

The calculation of Eq. requires the statistical average of products of four fermionic operators which is directly 
obtained using the Wick theorem and the expressions for the one-particle Green’s function that were derived in ||8l. 


<cJ(f')Q(0> = 2 r ^ME)%E(lt')^‘aEii,t) 

Qt 

<c,(f)c](0> = 2 f ^(^-f-(E))'i'aEii,t)^‘:Eij,t'). 

n 


One obtains. 


dEi f~* dEi ^ 


with the quantity l^EE'it) closely related to the initial current operator 


l,EE'it) = ^ 


(9) 

( 10 ) 

( 11 ) 


( 12 ) 


(13) 


Equations ( [T^ relates the typical output of a time-dependent simulation (right-hand side) to the correlation properties 
(left-hand side). To proceed, we focus on the total number of transmitted particles over a duration A, 

^A/2 

= I dtlit). (14) 

J-hji 

With these notations, the average and variance of are given by. 


<«/i) = 



faiE)NEE 


(15) 


var(n^) 




(16) 


with 


Nee' 



dt l,EE'it). 


(17) 


Eor a finite A, the above equations are well defined. However, as we shall see, the noise is totally dominated by the 
equilibrium noise which diverges at large A so that the numerical computation of the excess noise (the quantity which 
is usually measured) can become problematic. Ultimately, we would like to consider the limit where A is infinite. To 
proceed, we separate the equilibrium physics from the time-dependent one and introduce '^oEii, 0 which measures 
how the wave-function deviates from its stationary solution; 


'VaEii, t) = 'VUe t)e-'^' + t) 


(18) 


'kff£(h0 is actually the direct output of the techniques discussed in ijS). With these notations, one can perform the 
integration over time (in the limit of large A) and write the noise in terms of well behaved (converging) integrals. The 
variance of now reads 

var(«^) = cr^, A -H la-^ix + + Oil /A) (19) 
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where the three numbers cr^,, CTmix and cr^ are defined as. 


/.(£)(! - MEW^MOt ( 20 ) 

«',yS 

CTmix j ^ /.(£)(! -//3(£))Re[^££4,££(0)] (21) 

9 P dE dE' - 9 

tr' = 2 j I ^ “ //t(£')]|A^££f (22) 

in term of 

X OO 

dt (f) - 4,££-(0)e’''^-^')']. (23) 

oo 

The convergence of the time integrals is assured if the time-dependent perturbation is of finite duration. The first term 
of Eq. ( |T^ is simply the stationary noise of the problem (including the Johnson-Nyquist noise and some shot noise 
contributions) while the two last terms form the extra noise due to the time-dependent perturbation. Eq. is very 
convenient for practical calculations as the diverging stationary noise can be directly subtracted from the calculation 
and one does not suffer from numerical inaccuracies. Eq. ([T^ is the main result of this article. The two quantities 



Figure 1: Comparison of the amplitudes of the bare (black full line) with Yaji l^£,oP where the diverging term has been subtracted 

(red dash line) as a function of energy. 


that need to be integrated, Nee' (bare) and Nee' (after subtraction of the stationary state), are plotted for illustration in 
Fig.[T]for a voltage pulse sent through a simple one-dimensional infinite wire with a barrier in the middle (see 191 for 
the precise microscopic model). One observes that the interesting contribution of Nee' is completely overshadowed by 
the stationary contribution (that contains a term of the form sin(A£’)/£’ which eventually converges to a dirac function) 
so that practical calculations are greatly facilitated by the use of Eq. ( [T9| ) instead of Eq. ( [T6| . 

3.2. Application to a one-dimensional conductor 

The study of the noise associated to a voltage pulse in a quantum conductor was pioneered some time ago by 
the work of Levitov and coworkers ifTJl fT4ll but the first experiments were only performed very recently ifTSl . Let 
us quickly revisit this issue. We consider a one-dimensional wire with a central barrier of transmission (reflection) 
amplitude d (r), and we apply a voltage pulse V(t) on the left electrode while the right one is grounded, see Fig.j^for a 
schematic. Assuming a linear spectrum E - vk and neglecting the energy dependence of the barrier, the wave-function 
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Figure 2: Schematic of our one-dimensional wire connected to two electrodes. The gate voltage Vg creates a barrier with transmission (reflection) 
amplitude d (r). A voltage pulse is applied on the left electrode while the right one is grounded. 


on the right of the barrier for electrons coming from the left and right electrodes takes the form 

'f'i £(n, f) = (24) 

^/v 

„~ikn~iEt „ 

0 = -^ (25) 

Vv yv 

where the two indices L, R label the electrode (Left or Right), v - dEIdk is the group velocity and 0(f) = du eV{u)lh. 
Inserting these wave-functions into Eq. ( [T^ and defining e - E - E', we obtain 

5««(f,f') = ^/(£)(!-/(£-e))( 2 D 2 + 2 D(l-D)cos[ 0 (f)- 0 (f')])e'^''"''> (26) 

where D - \d'^ is the transmission probability of the tunneling barrier. Performing the successive integrals over the 
times f and f', and over the energy E one obtains for the variance, 

kgTeD^ r de 9 / f \ 

var(n^) = —-—A + D(1 - D) j (27) 

where K(e) = f dt jjjg amplitude of probability that the voltage pulse leads to a change e of the energy, ks 

is the Boltzmann constant and Tg the temperature. The second term of Eq. ( p7] ) is equal to Eq. (17) of Ref ifTSl up to 
the barrier factor D(1 - D). The partitioning of Eq. ( [27) i is however not very transparent as the (diverging) equilibrium 
noise is actually contained in both terms of the equation. In fact, in the absence of pulse and for infinite A, K{e) — 6(e) 
so that the second term of Eq. p7| contains the square of a dirac function and is ill defined. Let us now compute the 
three contributions cr^,, o-^ix and separately. The deviation from equilibrium of the wave-function takes the form 

'El E(n, t) = - l)e'^"-'^' (28) 

Vv 

'P«,£(«,f) = 0. (29) 


and one obtain after some algebra. 


cr„ = 


kgT.Dil-D) 


/-A/2 

' 

J-A/2 


dt cos(0(O) - A 


(T + 2,0~mix — ^(1 




kgTeDil-D) 


(30) 

(31) 

(32) 


One can easily check that the summation of all these terms gives back the result of Eq. ( [27| as it should. The advantage 
of this form is that the different integrals are now naturally regularized while the Johnson-Nyquist contribution is well 
separated from the rest. Equation ( |27) l was further analyzed by Levitov et al. in Ref. in and one obtains a variance 
of the form (at zero temperature). 


var(n^) - Di\ - D)n + D(l - D)a sin^( 7 rn) 


(33) 
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Figure 3: Voltage pulse across a simple banier. Variance of the number of transmitted particles at zero temperature as a function of the number of 
injected particles. The symbols correspond to the numerical data obtained by direct integration of the wave-function. The full line corresponds to 
Eq.i|33|with a ^ 3. The bamer transmission is set to D = 0.5. 


where a is a numerical factor and n - J dt eV{t)lh is the number of electrons sent by the voltage pulse. The oscillating 
behavior in sin^(7rh) has been extensively discussed in the context of recent experiments M- Similar oscillations are 
also predicted for the current itself in the context of interferometers i). 

To illustrate our numerical procedure, we calculate the variance of the number of particles transmitted in the case 
of Lorentzian voltage pulse, V(t) - Vp/(1 + {ItjTpf'), with Vp the amplitude and Tp the duration of the pulse (and 
h cx VpiTp). Fig.|^shows var(hp) as a function of h at zero temperature and for a semi-transmitting barrier D - 0.5. 
We get a perfect match with the analytical prediction. Fig.|^shows the noise var(hp) at various temperatures after 
subtraction of the equilibrium noise. For h = 1 the noise is minimum, i.e. equal to the shot noise created by an 
equivalent d.c. current. This is a peculiar feature of the Lorentzian shape which sends a unique electron (without 
creating additional electron-hole pairs that would cause extra noise). The amplitude of the oscillations shows a rapid 
decrease with temperature as shown in inset. 



Figure 4: Variance of the number of transmitted particles where the Johnson-Nyquist noise has been subtracted as a function of the number of 
injected particles. Avar is the difference of the noise for fi = 0.5 and h = \. Inset: Avar as a function of temperature in the unit of h/Tp. The dots 
are numerical data, and the continuous line is a fit obtained with Avar = a exp{~bhl(TphsTe)). The barrier transmission is set to D = 0.5. 


4. Conclusion 

We have presented an extension of the numerical technique of IS] that allows one to compute the quantum noise 
properties of a conductor. The method naturally separates the stationary noise from the extra contributions coming 
from time-dependent perturbations such as voltage pulses. Once the wave-functions are known for different values 
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of the energy, the calculation of the noise simply amounts to a numerical integral which can be performed as a 
post process. As an application we have considered the simple case of a one-dimensional system but the numerical 
technique is in no way limited to this case. In particular, the large class of systems which have already been considered 
for the time-dependent simulations (which includes the quantum Hall effect, electronic interferometers, conventional 
and topological superconductors...) can now be studied for their noise properties. Higher moments, correlations or 
fluctuations of other observables (such as charge fluctuations) can also be obtained following the same lines. 
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